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D . Abstract 



p. 

The lightest neutraUno is a viable dark matter (DM) candidate. In this pa- 
^ ■ per we study indirect detection of the wino-like neutralino DM using positrons 

' and antiprotons from the annihilation in the galactic halo. When the mass is 

around 2 TeV, which is favored from the thermal relic abundance, the non- 
perturbation effect significantly enhances the annihilation cross sections into 
positrons and antiprotons. We find that the positron and antiproton fluxes 
with energies larger than 100 GeV may become larger than the expected back- 
grounds. Since the positron fiux is less sensitive to the astrophysical param- 
eters, the detection may be promising in the upcoming experiments such as 
PAMELA and AMS-02. We also find the wino-like neutralino DM with mass 
around 2 TeV is compatible with the HEAT anomaly. 



1 Introduction 



The existence of the cold dark matter (CDM) has been confirmed by the WMAP 
measurement of the cosmic microwave background [l]; ncDuh^ = O.lUtoms 13 • 
However, the nature of the dark matter (DM) still remains a mystery. Weakly 
Interacting Massive Particles (WIMPs) are viable candidates for the DM since their 
thermal relic abundances are naturally within the observed range 

A well-studied representative of WIMPs is the lightest neutralino in supersym- 
metric extensions of the standard model 0]. Neutralinos are composed of bino, neu- 
tral wino and neutral Higgsinos, which are superpartners of the U(l)y and SI]{2)l 
gauge bosons and the Higgs bosons, respectively. In most supersymmetric models, 
the lightest neutralino is the lightest supersymmetric particle (LSP) and is stable 
due to the R-parity conservation. The constituent of the neutralino depends on su- 
persymmetry breaking models. For example, the neutralino is bino-like in a wide 
region of the parameter region of the minimal supergravity model. In the anomaly 
mediated supersymmetry breaking model (AMSB) j^, the neutralino is wino-like 
because gaugino masses are proportional to beta functions of the gauge coupling 
constants. The wino-like neutralino has larger coupling than the bino-like one, so 
that the wino-like neutralino DM has larger prospects for detection. 

Various experiments have been performed or are planed in order to detect the 
neutralino DM directly or indirectly. The direct detections are to measure the recoil 
energy which the neutralino may deposit as it crosses a terrestrial detector |H] . The 
detection rate depends on the cross section for the elastic-scattering of the neutralino 
with target nuclei. On the other hand, the indirect ones are to detect the anomalous 
cosmic rays produced in the neutralino annihilation. Detectors are designed to ob- 
serve high-energy neutrinos from the earth or the sun, gamma-rays from the galactic 
center, and antimatter cosmic rays from the galactic halo, which are generated from 
the neutralino annihilation. 

In this paper, we consider indirect detection of the wino-like neutralino DM by 
positrons and antiprotons in cosmic rays. It is pointed out in Refs. j71IH] that the cross 
sections for the wino-like neutralino annihilation into gauge bosons are enhanced 
compared with those at the tree-level approximation when the mass is larger than 
about 1 TeV. This is due to a non-perturbative effect by the electroweak interaction, 
which appears in a non-relativistic limit of the wino-like neutralinos. Especially, 
when the mass is around 2 TeV, which is favored from the thermal relic abundance of 
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the wino-like neutralino [Oj, the annihilation cross sections are significantly enhanced 
by the resonance effect. 

The enhancement of the annihilation cross sections raises the possibilities of the 
indirect detection of the wino-like neutralino DM. In Ref. 0, the gamma-ray flux 
produced by the wino-like neutralino annihilation in the galactic center is evaluated, 
and it is found that the the sensitivity for the wino-like neutralino DM is enhanced. 
In this paper, we evaluate the positron and antiproton fluxes from the wino-like neu- 
tralino annihilation in the galactic halo, including the non-perturbative effect. We 
flnd that, for the neutralino with mass around 2 TeV, the positron and antiproton 
signals also exceed the backgrounds. These fluxes will be measured with unprece- 
dented accuracies by the upcoming experiments such as PAMELA ^IflJ and AMS-02 
jTTj. Especially, the measurement of the positron flux may be more promising for 
detection of the wino-like neutralino with mass around 2 TeV, since the predicted 
positron flux is less sensitive to the astrophysical parameters responsible to the prop- 
agation or the DM halo proflle. 

This paper is organized as follows. In section 2, we review the non-perturbative ef- 
fect on the wino-like neutralino annihilation cross sections. In section 3, the positron 
flux from the annihilation in the galactic halo is evaluated using the diffusion model. 
Here, we compare the predicted signal positron flux with the expected background, 
and discuss the sensitivities of the future experiments to the heavy wino-like neu- 
tralino DM. The HEAT anomaly [12J is also discussed. In section 4, we investigate the 
antiproton flux from the wino-like neutralino annihilation. The expected background 
and the future prospect are also discussed. Section 5 is devoted to conclusions. 



2 Non-perturbative effect on wino-like neutralino 
annihilation 

The wino-like neutralinos annihilate mainly into W bosons due to the SU(2)i gauge 
interaction. The annihilation process is mediated by t-channel wino-like chargino 
exchange at tree-level, and the cross section is given by 

27ra2 

= — r ' (1) 

where v is the relative velocity of the neutralinos, a2 is the 811(2)^ gauge coupling 
constant and m is the wino mass. 
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We take a non-relativistic limit {v <^ 1) in Eq. however, the tree-level approx- 
imation in the limit is not valid for the wino-like neutralino heavier than ~ mvi//a2- 
Here, my/ is the W boson mass. This is due to the threshold singularity caused 
by the mass degeneracy between the wino-like neutralino and chargino, and the 
higher-order contributions should be included in the case. The dominant contribu- 
tion to the scattering amplitude at 0(0:2"^^'') comes from ladder diagrams, in which 
n gauge bosons are exchanged. When the mass difference between the wino-like 
neutralino and chargino is negligible, the n-th ladder diagram is suppressed by only 
{a2'm / mw)'^"^^ compared with the leading-order one ^Hj- Thus, when m is larger 
than ~ mwloi2, we need to resum the diagrams at all orders. In other words, we 
have to include the non-perturbative effect for obtaining the reliable annihilation 
cross section. 

The resummation of the ladder diagrams has a following interpretation. Since the 
wino mass is much heavier than that of the W boson, the wino-like neutralinos feel 
the long-range force induced from the W boson exchange. Due to the force, the wave 
function of the neutralino pair is significantly modified from the plane wave before 
the annihilation into W^W~ bosons. As shown in Ref. 0, the bound states, which 
are composed of the neutralino and chargino pairs, appear due to the long-range 
force if the wino mass is large enough. Especially interesting, a bound state has the 
binding energy almost zero when the wino mass is close to ~ 2, 8, • • • TeV. In those 
cases, the wino-like neutralino annihilation cross section in a non-relativistic limit 
is enhanced by several orders of magnitude compared to that of the tree-level cross 
section due to the resonance. 

In Fig. the annihilation cross sections into W^W~ and ZZ bosons are shown 
as functions of the wino mass. These figures are plotted using fitting formulae for 
the wino-like neutralino annihilation cross sections given in Ref. |H]. When the 
mass difference between the wino-like neutralino and chargino is much smaller than 
Oi2mw, which is a typical potential energy due to the electroweak interaction, the 
cross sections are less sensitive to the value of the mass difference. In this paper, the 
mass difference is set to be 0.1 GeV for definiteness. For heavy wino-like neutralino, 
this mass difference is dominated by the radiative correction, and it is 0.1 — 0.2 GeV 
in most of the parameters region [8 . This is because the tree-level contribution to the 
mass difference is suppressed by (mv^^/MsusY)"^ unless the wino mass is accidentally 
finely tuned to the Higgsino mass. The mixing between the wino and Higgsino 
components is also suppressed by (mviz/Msusv)- Thus, we ignore the mixing in the 



3 



Annihilation cross section into WW 




m ( TeV ) 



Annihilation cross section into ZZ 



10 



-23 



10 



-25 



10 



-27 



0.2 



1 

m f TeV ) 



(cm-'s 


j 

u 


6m = 0.1 GeV 





10 



Figure 1: Cross sections, av, of the annihilation of the wino-like neutralinos into 
W^W~ (left figure) and ZZ (right figure) in a non-relativistic limit. The mass 
difference between the wino-like neutralino and chargino is set to be 0.1 GeV. For 
comparison, the cross sections at the leading order in perturbation are shown as 
dashed lines. The bound state resonances appear around 2 TeV and 8 TeV. 



following. 

As shown in the figure, the annihilation cross section into ZZ is also enhanced for 
m ^ 1 TeV in addition to that into W^W~ , and it becomes comparable to that into 
W^W~ . The cross sections into 77 and 7Z also have a behavior similar to that into 
ZZ. The annihilation channels into ZZ, 7Z and 77 come from one-loop diagrams 
in the perturbation, and the cross sections are suppressed. However, the transition 
between the neutralino pair state and the chargino pair state is not suppressed 
due to the non-perturbative effect for m ^ 1 TeV, so that the cross sections are 
enhanced. When evaluating the positron and antiproton fluxes from the wino-like 
neutralino annihilation in the galactic halo, we need to include the contribution of 
the annihilation into Z bosons, in addition to that into W bosons. 

If the relic abundance of the wino-like neutralino in the universe is explained by 
the thermal scenario, the mass consistent with the WMAP observation is around 
2 TeV [9J. It is intriguing that this value is coincident with the mass corresponding 
to the resonant annihilation as shown in Fig. ^ On the other hand, the wino-like 
neutralino DM is also produced by non-thermal processes such as the moduli decay 
jl4H 15j. Furthermore, the late time entropy production by, for example, the thermal 
inflation ^S] may decrease the amount of the DM. In these cases, the mass of the 
wino-like neutralino consistent with the DM observations may be deviated from 
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2 TeV. 

In this paper, while the heavy wino-hke neutrahno with mass around 2 TeV 
is noticed, we discuss the positron and antiproton signatures from the neutrahno 
annihilation without peculiar masses specified for completeness. Thus, we assume 
that the wino-like neutrahno is dominant constituent of the CDM in the present 
universe, and exists in the halo of our galaxy with appropriate mass density in the 
following. 

3 Positron signature of wino-like neutralino dark 
matter 

In this section, we evaluate the positron flux from the wino-like neutralino annihi- 
lation in the galactic halo. In the evaluation of the signal flux in the vicinity of the 
solar system, we need to consider the propagation of positrons through the galaxy, 
in addition to the production rate of the positrons from the annihilation in the halo. 
We discuss these in order, and show the sensitivities of the upcoming experiments, 
such as PAMELA and AMS-02, to the positron signal by comparing the expected 
background originated from the secondary production of the cosmic rays. The HEAT 
anomaly is also discussed. 

3.1 Production rate of positrons from dark matter annihi- 
lation 

The production rate of positrons from the neutralino DM annihilation in the galactic 
halo is given as 

g(£;,r) = in^(r)5]M,(^)^ , (2) 

where n is the number density of the neutralinos in the galactic halo, {av) f is 
the annihilation cross section into the final state /. The fragmentation function 
{dNe+ 1 dE) f represents the number of positrons with energy E, which are produced 
from the final state /. The coefficient 1/2 comes from the pair annihilation of the 
identical particles. 

As discussed in the previous section, the wino-like neutralinos annihilate into W 
and Z bosons. Positrons are produced through the leptonic and hadronic cascade 
decays of the weak gauge bosons, for example, W'^ — > e'^u, — > ^i^v — > e^vvv 
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Figure 2: Fitting functions of the fragmentation functions {dNg+ /dx)ww ^-nd 
{dNf.+ /dx)zz (solid lines) and HERWIG Monte-Carlo results in cases of m = 0.2, 0.5, 1, 1.5, 
and 2 TeV. 



or — >• hadrons —> tt /i"*" —> e"*". These cascade decay processes for producing 
positrons are encoded into the fragmentation functions {dNe+ / dE) f (/ = WW and 
ZZ). We ignore the contribution from the annihilation into Z7, since the contri- 
bution is less than about 10%. We evaluate the fragmentation functions using the 
HERWIG Monte-Carlo code [TT and derive the fitting functions as follows, 

= exp [WW{\n{x))] , (^) = exp [ZZ(ln(x))] , (3) 
\ ax J T^y^ \ ax J ZZ 

where x = E/m and the functions, PFPF(x) and ZZ{x), are given by 

WW{x) = -2.28838 - 0.605364a; - 0.287614x2 - 0.762714x^ 

-0.319561s^ - 0.0583274a;^ - 0.005035552;'^ - 0.00016691x^ , 
ZZ{x) = -2.75588 - 0.45725x - 0.141373^^ - 0.905392x^ 

-0.444098x^ - 0.0936451x^ - 0.00942148a;'^ - 0.000369777a;^ . (4) 

In Fig. 121 the fragmentation functions from the HERWIG code and the fitting 
functions are depicted. The results of Monte-Carlo simulations are shown for cases 
of m = 0.2, 0.5, 1, 1.5, and 2 TeV. The fitting functions are shown as solid lines and 
agree well with the simulation data with the range m > 300 GeV and x > 10"'^. It 
is found that the slopes of the fragmentation functions are changed around x ~ 0.2. 
The positrons with lower energy {x < 0.2) comes from the hadronic cascade decay 
process ^HJ; while those with higher energy (x > 0.2) are produced more directly 
from the leptonic weak boson decays. 
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Next, we discuss the DM number density in the galactic halo. The number density 
is derived from the DM halo mass profile p(r') through the equation n{f) = p{f)/m. 
The halo mass profile is determined by observations of the rotational velocity of the 
galaxy and the motions of the dwarf galaxies with help of the A^-bodies simulations, 
while several models for the DM profile are proposed. In this paper, we use the 
isothermal halo model, which is given as 

where r = |r| is the distance from the galactic center, 0.43 GeV/cm^ is the local 
halo density (the mass density in the vicinity of the solar system), 2.8 kpc is the core 
radius of the galaxy, and 8.5 kpc is the distance between the galactic center and the 
solar system. 



3.2 Propagation of positrons in the galaxy 

Once positrons are produced by the DM annihilation, they travel in the galaxy under 
the influence of the tangled magnetic field. Since the typical strength of the magnetic 
field is a micro Gauss, the gyroradius of the positron is much less than the galactic 
radius. Thus, the propagation can be treated as a random walk, and only some 
portion of the positrons can reach to the earth. 

There are some models for the propagation. Among those, we use the 'diffusion 
model' in which the random walk is described by the diffusion equation, 

^/e+(i?,f) = K{E)V'fME,f) + A [b{E)U{E,f)] + Q{E,r) , (6) 

where fe+{E^r) is the number density of positrons per unit energy, E is the energy 
of positron, K{E) is the diffusion constant, b{E) is the energy loss rate, and Q{E, r) 
is the source (positron injection) term discussed in the previous section. The flux 
of positrons with high energy {E ^ me) in the vicinity of the solar system is given 
from fe+{E, f) as 

$e+(i?) = ^/e+(i?,ro) , (7) 

where c is the velocity of light and Tq represents the coordinate of the solar system. 

The diffusion constant K{E) in Eq. ® is obtained by the simulation of cosmic 
rays, in which the diffusion model is used. In particular, the Boron to Carbon ratio 
B/C is an important quantity for the simulation. By comparing the measurement 
of B/C in the cosmic rays and the result of the simulation, the diffusion constant is 
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evaluated. For the calculation of the positron flux, we use the value in Refs. jTHl ED] , 

K{E) = 3 X 10^7 [3°-^ + {E/1 GeV)°-^] (cm^s-^) , (8) 

where the form of K{E) affects low energy positron flux, while high energy one which 
we are interested in is almost independent of the choice of this parameter. 

The positrons lose their energies by the inverse Compton scattering with cosmic 
microwave radiation (and infrared photons from stars) and the synchrotron radiation 
with the magnetic field during the propagation in the galaxy. Therefore, the energy 
loss rate b{E) is determined by the photon density, the strength of the magnetic field 
and the Thomson scattering cross section. We use the value of b{E) in Refs. [ini EI] , 

b{E) = 10-^%E/1 GeVf (GeVs-^) . (9) 

It is plausible that the positrons from the DM annihilation are in the equilibrium 
in the present universe, and hence the number density /e+ {E, r) is obtained by solving 
Eq. © with the steady state condition dfe+/dt = 0. Furthermore, we impose the 
free escape boundary condition, namely the positron density drops to zero on the 
surface of the diffusion zone. The positrons coming from the outside of the diffusion 
zone are negligible, and the positrons produced inside the diffusion zone contribute 
to the flux around the solar system, since they are trapped due to the magnetic field 

m 

It is usually assumed that the diffusion zone is a cylinder and that its half-height 
and radius are L ~ (2 — 15) kpc and R = 20 kpc, respectively. We fix L = 4 kpc 
in the evaluation of the positron flux. However, high-energy positrons, which we 
interest, only come from within a few kpc of the solar system as will discussed later. 
Hence, the positron flux is weakly dependent on the choice of the parameters of 
the diffusion cylinder. A detailed method for solving the diffusion equation (jH)) is 
presented in Appendix A. 

Here we discuss the effect of the solar modulation on the positron flux. The flux 
given by Eq. (|7j) is not exactly one to be measured on the top of atmosphere. The 
spectrum of the interstellar flux in Eq. ((7j) is modified due to interaction with the 
solar wind and the magnetosphere. However, the effect is not so important when the 
energy of the positron is above 10 GeV. Furthermore, the solar modulation effect is 
removed in the positron fraction, that is a ratio of positron to the sum of positron 
and electron fluxes, e"'^/(e+ + e^). Thus, we present our result mainly in terms of 
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the positron fraction. 

3.3 Background fluxes of positrons and electrons 

Positrons in the galaxy are injected by not only the DM annihilation but also the 
scattering of cosmic-ray protons with the intersteller medium. (See e.g. 23^^.) The 
flux of these positrons is calculated by simulations, in which the diffusion model is 
also used. The results agree with the measurements of the low-energy positron flux 
in the cosmic rays [23] • 

Since we can not distinguish the signal positrons, which originate from the DM 
annihilation, from those background positrons in measurements, we need to know 
the background positron flux. The background electron flux is also required for 
predicting the signals in terms of the positron fraction. In this paper, we use the 
fitting functions of these background fluxes, which are obtained by the cosmic ray 
simulations [T^ . 

= (GeV-^cm'2 -1 -1) 

^(sec).^s O^TOE^ fGeV-^cm-^s-^sr-i) 

1 + 110^1-5 + 600^2.9 + 580^4.2 l^^^ cm s sr J, 

^^T\e) = (GeV-^cm-^s-^sr-i) , (10) 

^+ ^ ^ 1 + 650^2.3 + 1500^4-2 ^ ^ ' ^ ^ 

where E is in unit of GeV. The first one, is the flux of the primary electrons. 

These electrons are considered to be produced by the shock wave acceleration in 
supernovae. On the other hand, the second and third ones, $^1''^ and (!>^T\ are 
the secondary electron and positron fluxes, respectively, which are produced by the 
collisions of cosmic ray protons and helium nuclei with hydrogen and helium of 
interstellar medium. 

3.4 Positron signature from dark matter annihilation 

In this section, we present the signature of the positrons from the wino-like neutralino 
DM annihilation. The positron flux from heavy DM annihilation (m > 1 TeV) is 
usually expected to be small. This is because the source injection Q scales as oc 
due to the mass dependence of the cross section (oc m~^) and that of the number 
density squared (oc m~2). However, the mass dependence of the cross section is very 
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Figure 3: (Interstellar) positron flux from the wino-like neutralino annihilation. The signal 
fluxes for the wino mass m = 0.3, 0.6, 1, 1.5, 2, 2.5, and 3 TeV are shown as solid lines. 
The expected background flux of positrons from the cosmic ray simulation is also shown 
as a dotted line. 



different from the ordinary one when the DM is the wino-like neutralino as discussed 
in the previous section. Furthermore, the cross section is enhanced by several orders 
of magnitude when the neutralino has the mass around 2 TeV. Thus, the positron 
flux is expected to be large in this case. 

First, we show the positron flux from the wino-like neutralino annihilation in 
Fig. 121 In this figure, the signal flux is shown as solid lines. The wino mass is taken to 
be m = 0.3, 0.6, 1, 1.5, 2, 2.5, and 3 TeV. For comparison, the expected background 
flux of positrons from the cosmic ray simulation is also shown as a dotted line. The 
effect of the solar modulation is not included, and thus the spectrums below 10 GeV 
have uncertainties. However, since the high-energy positron spectrum is important 
for the discrimination of the signal from the background as indicated in Fig. El the 
uncertainties from the solar modulation are not serious. 

When the wino mass is around 300 GeV, the signal flux is comparable to the 
background flux in the energy range 100 GeV<E< 300 GeV. Furthermore, the 
signal flux for the mass around 2 TeV also exceeds the background one in the energy 
range E > 100 GeV. The latter comes from the resonant DM annihilation. It is also 
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Figure 4: (Left figure) Positron fraction, 6^/(6"*" + e~), as a function of positron energy 
E in the wino-like neutralino DM. For comparison, the expected background positron 
fraction, the positron data HEAT 94-95 and HEAT 2000 are also shown in this figure. 
(Right figure) Contour plot of the ratio between the positron fractions including positrons 
from the DM annihilation and without it (that is the background positron fraction) in a 
{E, m) plane. 

noticed that a bump appears in each signal spectrum at around m/2. The positrons 
with energy above the bump come from the direct decay of weak gauge bosons, while 
those with energy below the bump are produced mainly by the hadronic cascade 
decay of the gauge bosons. 

Next, we consider the positron fraction calculated from the positron flux in the 
Fig. island the expected background ones in Eqs. (fTUj) . The result is shown in Fig. HI 
In the left figure, the positron fraction is depicted as a function of positron energy 
for several wino masses. The choice of the mass is the same as that in the Fig. El 
The expected background positron fraction, the positron data HEAT 94-95 ^21 and 
HEAT 2000 [21] are also shown in this figure. In the right figure of Fig. |31 the ratio 
of the fraction including positrons from the DM annihilation to the background one 
is depicted as a contour plot in a {E, m) plane. From these figures, it is clear that the 
signature becomes more significant for high-energy positrons. In particular, there is 
a large difference between the expected signal and the background when the wino 
mass is a few hundred GeV or around 2 TeV. 

Here, we address the HEAT experiment ^3 121]; which reported the positron 
excess from the expected background. The spectrum of the observed fraction is 
almost flat around 0.06 in a energy range 4 GeV ^ -E < 20 GeV. The positron 
fractions for both m = 300 GeV and 2 TeV in the figure are consistent with it within 
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Figure 5: Sensitivities of the upcoming experiments. The positron fraction, e^/ (e"*" + e~), 
for m = 2 TeV and that of the background are shown as solid and dotted lines, respectively. 
The error bars in the figure correspond to the statistical errors projected for the PAMELA 
and AMS-02 experiments after three years of observations. 

the experimental error. 

In addition, the effect of the inhomogeneity in the local DM distribution on the 
positron flux is recently discussed, whose existence is supported by the A^-bodies 
simulations. In these arguments, the positron flux from the DM annihilation is 
enhanced if there are clumps of the DM in the vicinity of the solar system. The 
effect is parametrized as a boost factor {BF) 25j, which is defined by a ratio of 
the signal fluxes with inhomogeneity and without inhomogeneity. The boost factor 
may reach ~ 5 when the inhomogeneity exists, while the factor is equal to one if 
the DM is distributed homogeneously. Thus, the wino-like neutralino with the mass 
~ 300 GeV or 2 TeV can explain the HEAT result quite naturally. It is amazing 
that the wino-like neutralino with 2 TeV naturally accounts for not only the DM 
abundance thermally but also the HEAT anomaly. 

Next, we discuss the potential of the upcoming PAMELA TIT and AMS-02 [TT] 
experiments, which have good sensitivities in a broad region of positron energy 
10 GeV ^ E < 270 GeV, might detect the signal from the wino-like neutralino 
DM annihilation. We estimate the sensitivities of those experiments following the 
method in Ref. ^7\. In Fig. we show the sensitivities of the PAMELA and AMS- 
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02. The positron fraction, /{e^ + e~), for m = 2 TeV and that of the background 
are shown as sohd and dotted hnes, respectively. The error bars in the figure corre- 
spond to the statistical errors projected for the PAMELA and AMS-02 experiments 
after three years of observations. As shown in this figure, positrons with energy of 
some tens of GeV will be clearly discriminated from the background. 

Finally, we discuss other uncertainties of the signal flux. First, in the case of 
the positron propagation with high energy, we do not have to worry about un- 
certainties from the thickness of the tangled magnetic fleld (L). This is because 
high-energy positrons we observe are produced within a few kpc around the solar 
system. Positrons far from the earth lose their energies during the propagation, and 
consequently they contribute to the low-energy part of the flux. The distance in 
which positrons travel without signiflcant energy loss is typically 



Thus, the positron flux at high energy does not suffer from the uncertainties of the 
thickness L (because L > a few kpc). 

Second is the DM distribution in the halo. We have assumed the isothermal 
halo in Eq. ^ in the above. Various DM halo models are proposed from the A^- 
bodies simulations, however, the high-energy positron flux from the DM annihilation 
is considered to be almost independent of the choice of the halo model. The main 
difference among the halo models appears in the galactic center. However, the high- 
energy positrons produced around the galactic center can not reach to the earth, and 
the positron flux has little ambiguity from it around the solar system. 

4 Antiproton signature from wino-like dark mat- 
ter annihilation 

The antiproton flux from the wino-like neutralino DM annihilation is discussed in 
this section. The method for calculation of the flux is essentially the same as that of 
the positron flux. First, the antiproton injection in the galactic halo (source term) 
and the propagation of the antiprotons are discussed, and the antiproton flux from 
the wino-like neutralino annihilation in the galactic halo is evaluated. The antiproton 
background originated from the cosmic rays is also discussed. 
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panel) . 





i = i 


J = 2 


J = 3 


i = 4 


i = 1 


480.0 


0.26 


9.6x10-"' 


2.27 


i = 2 


2.17 


0.05 








i = 3 


-8.5 


-0.31 








j = 4 


-0.33 


-0.075 


-1.5 X lO-'' 


0.71 



W process (left panel) and ZZ one (right 



4.1 Production rate of antiprotons from dark matter anni- 
hilation 

Antiprotons from the wino-like neutralino annihilation are also produced through 
the cascade decay of weak gauge bosons. The difference between the antiproton 
and the positron production rates (source terms) appears only in the fragmentation 
functions, and then the antiproton production rate is given as 

Q(T,r) = in^(f)5^M,(^) . (12) 

where T(= E — m) is the kinetic energy of antiproton and {dNp/dT) f is the frag- 
mentation function. 

As in the case of the fragmentation functions for positrons, the functions for 
antiprotons are calculated from the Monte-Carlo simulation. In this paper, we use 
the simple parametrization in Ref. |2H1 which fits the result of the PYTHIA Monte- 
Carlo code [21, 

^) =iPixP'+P2\log,,xn~' , (13) 

where x = T/m. The parameters pi in the above equation depend on the neutralino 
mass in addition to the annihilation channels, and they are given as 

Pi{m) = {annf^^ + aism"^"^)"^ . (14) 

The values of the coefficients, listed in Tabled for the W~^W process (left 

panel) and the ZZ one (right panel). The parametrization for the fragmentation 
functions is valid for the neutralino mass in the range (50 - 5000) GeV. We dropped 
quark processes such as ti and hh since the annihilation cross sections are very small 
due to heavy squark masses and helicity suppression. 
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4.2 Propagation of antiprotons in the galaxy 

In order to treat the propagation of antiprotons, we use the diffusion model as in the 
case of positrons. The diffusion equation describing the propagation is written as 

KpiT)V'fpiT,r) - ^ iVciz)fp{T,r)) - 2hS{z)T,,Jp{T,r) 

+Q{T,r) + Q"^^\T,f) = 0, (15) 

where fp{T, r) is the number density of antiprotons per unit energy. The steady 
state condition {dfp/dt = 0) is assumed as discussed in the positron case. For 
the evaluation of the equation, we use the cylinder coordinate. The interaction of 
antiproton with matter is confined on the galactic plane, which is expressed as the 
infinitely thin disk with radius R = 20 kpc at the z = 0. The diffusive halo is the 
cylinder with radius R = 20 kpc and the half-height L. The boundary condition for 
solving the equation is taken to be the same as the positron case. 

The diffusion equation (|T3|l is essentially the same as that in Eq. ©• However, 
there are some differences, for example, the energy-loss term does not appear in 
Eq. ()15|). This is because protons are much heavier than electrons, so that we can 
neglect the energy-loss due to the scattering with background photons. The other 
differences are the term related to the convective wind (second term), the interaction 
term with matter in the galactic plane (third term) and the tertiary antiproton term 
(last term). These three terms are not so important when we consider the antiproton 
flux with high energy (T > a few GeV). We include these terms in the diffusion 
equation for completeness. 

The diffusion coefficient Kp is determined by the Boron to Carbon ratio B/C in 
the cosmic rays, which is the same as the positron case. For the calculation of the 
antiproton flux, we parametrize the diffusion constant as Refs. \60\ I31j. 

Kp{n) = Kopn' , (16) 

where the diffusion coefficient is assumed to be constant within the diffusion zone. 
The variable TZ is called the rigidity, which is defined by the momentum of the 
particle per unit charge 71 = p/Z. For the values of 6 and Kq, we use the parameter 
sets in Table |21 These values are favored from the B/C analysis [22] • 

The second term in Eq. (fT^ . d(yc{z)fp)/dz, is not included in the equation for 
the positron fiux. This term is related to the convective wind, which represents the 
movement of medium responsible for the diffusion. The direction of the wind is 
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case 


6 


i^o(kpc7Myr) 


L (kpc) 


Vc{km/s) 


max 


0.46 


0.0765 


15 


5 


med 


0.70 


0.0112 


4 


12 


mill 


0.85 


0.0016 


1 


13.5 



Table 2: Astrophysical parameters compatible with the B/C analysis [HSj. Three cases 
give the maximal, median and minimal signal antiproton fluxes. 



assumed to be perpendicular to the disc plane, and the velocity Vc{z) is constant 
throughout the diffusive volume, 

Vc{z) = mz) - 1) K , (17) 

where the value for Vc is given in Table El 

Next one is the third term in Eq. (fT^ . —2h6{z)Ts,^^fp, which represents the 
annihilation between antiproton and interstellar proton in the galactic plane. The 
parameter h in the term is the half- height of disk and set to be 100 pc L), while 
Fann is the annihilation rate between antiproton and proton, 

Tann = (r^H + 42/3nHe)a|^"^p- , (18) 

where Vp is the velocity of antiproton, nn denotes the hydrogen number density (~ 
1 cm~'^), and nne is the helium number density which we assume to be 7% of rin 
IHHl . The factor 0^ arises from a geometrical approximation [31] . The annihilation 
cross section between antiproton and proton, Cpp^, is given by 



661 (1 + 0.0115 r-°-^^^ - 0.948 T'^ o^^^) mb , T < 15.5 GeV 
■PP ^ 1 36 r-°-5 mb , T> 15.5 GeV 



where T is in unit of GeV. This interaction dominates over inelastic interactions 
at low energy. Hence, the flux of antiprotons with low energy is decreased by the 
annihilation. 

For higher energy antiprotons (T > 10 GeV), the inelastic interaction is not 
dominated by annihilation, however, the non-annihilating scattering is important. 
The interaction lowers energies of antiprotons, T' to T(< T'). These antiprotons are 
called tertiary antiprotons. We include this effect in Q^^^^{T,r), which is given by 



gt-t(T,f) = (nH + 42/3nHe) 

cm non— ann/T-i/ 
Opp yi 



X 



^ ' v'pfpiT, r)dr - a^^-^^^^{T)vpfp{T, r) 



.(20) 
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The first term in the bracket is the contribution to the antiproton flux with energy 
T from the inelastic scattering of antiprotons with energy larger than T, while the 
second term compensates it so that the total antiproton number is not changed in 
this process. Here, ap°^'^'^^{T) is given as the difference between the total inelastic 
cross section o"?f^ and the annihilation cross section cj'^^. The total inelastic cross 
section is given in Ref. [^3] as 

af^\T) = 24.7 (1 + 0.584T-°-i^5 + 0.856T-^-^^^) (mb) , (21) 

where T is in unit of GeV. 

The number density of antiprotons fp is obtained by solving the diffusion equation 
(fT^ . We can solve this equation fuU-analj^ically [221 EOl- The detailed expression of 
the solution is presented in Appendix B. After solving the equation for the number 
density, the interstellar flux of antiprotons from the DM annihilation in the vicinity 
of the solar system is obtained as 

^,s = ^fp{T,fe)- (22) 

Here, we discuss the effect of the solar modulation on the antiproton flux. This is 
important for antiprotons with low kinetic energies (< 3 GeV). Using the force fleld 
approximation, the flux of antiprotons on the top of atmosphere $toa is obtained 
from the interstellar flux <l>is as 

where pTOAj^j^TOA-j ^^^^ p^^(T^^) are momentums (kinetic energies) of antiproton on 
the top of the atmosphere and in the interstellar, respectively. The solar modulation 
parameter varies according to the 11 years solar cycle. This parameter takes a 
value from about 500 MV at the minimum solar activity to 1.3 GV at the maximum 
solar activity. Larger lowers the antiprotons flux on the top of atmosphere flux at 
low energy. 



4.3 Background flux of antiprotons 

In this section, we discuss the antiproton background flux. The antiprotons are 
produced as the secondary products of cosmic rays by the nuclear reaction with the 
interstellar gas in the galactic disk. The main contribution to the antiproton flux 
comes from the collision between the cosmic ray protons and the interstellar hydrogen 
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gas. Again, the production phenomena are described by the diffusion equation. We 
solve the diffusion equation and calculate the background ffux. Since the concrete 
formalism for obtaining the flux is very complex, we mention only the strategy for the 
calculation here. The antiproton background is also discussed in Refs. fIE\ ^^ ^7\. 

While the interstellar primary proton flux is required to evaluate the background 
antiproton flux, it is impossible to measure it directly. However, it is obtained by 
solving the diffusion equation under an assumption of the source function. The 
primary protons are believed to be produced by supernovae. Hence, the proton 
source term with a few undetermined parameters is assumed, and the interstellar 
proton flux is obtained by solving the diffusion equation with this source term. In 
this case, the parameters in the source term are fixed by comparing the evaluated 
flux with the observed cosmic rays on the earth in the measurements such as BESS 
PH] and AMS [321 • The fitting function for the primary proton flux derived as above 
is given in Ref. [31]. We use it in our evaluation of the background antiproton flux. 

Next, the antiproton flux is evaluated from the primary proton flux by solving the 
corresponding diffusion equation. The equation is the same as that in Eq. (fT^ except 
for the source term. Since the antiprotons are produced by the nuclear reaction 
between the cosmic rays and interstellar gas, the source term is given by the proton 
flux and the cross sections for the reactions. 

The antiprotons are dominantly produced by the process p + H — p + X. In the 
rest frame of the hydrogen atom, the kinetic energy threshold for the incident proton 
to produce secondary antiprotons is GrUp. Furthermore, the number density of the in- 
cident proton decreases as energy increases. As a result, the spectrum of antiprotons 
from this process has a peak at a few GeV. In addition to this process, we include 
the inelastic collision between proton and helium, p + He — > p + X, for generating the 
secondary antiprotons. The process with the helium contributes to the antiproton 
flux sub-dominantly in the most energy range. However, the antiprotons from the 
process are a dominant component at low energy with the tertiary antiprotons (T < 
0.1 GeV). Thus, the source term for the secondary antiproton turns out to be 

POO 

g(T) = 2 / c/T'47r$J,P""^)(r') 

The factor 2 comes from the fact that the antiprotons are produced from the antineu- 
tron decay in addition to the direct production of antiprotons. The threshold energy 
Tth is Grrip and ^^^^^^ is the proton flux. The differential cross section da{T', T)/dT' 
is for the production of an antiproton with energy T from an incident proton of 



(icT. 



dV 



dV 



iT\T) 



■ (24) 
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Figure 6: Antiproton flux from the wino-Iike neutralino annitiilation on the top of 
atmosphere as a function of antiproton kinetic energy. The wino mass is taken to be 
m = 0.3, 0.6, 1, 1.5, 2, 2.5, and 3 TeV. The solar modulation parameter is set to 
be 500 MV. For comparison, the background flux is depicted as a dotted line. 

energy T'. The cross sections are given in Ref. |40j . 

Within uncertainties of the observations, the obtained flux for the antiproton 
background is consistent with the results by BESS HI], AMS ^ and CAPRICE P] . 
which observe the low-energy antiprotons ((0.2-50) GeV). We use the background 
flux for estimating the antiproton signature from the DM annihilation. 

4.4 Antiproton signature from dark matter annihilation 

Now we are in a position to discuss the antiproton signature from the wino-like 
neutralino annihilation in the galactic halo. We calculate the antiproton flux from 
the neutralino annihilation by solving Eq. (fT3j) . and compare the result with the 
background flux discussed in the previous section. 

In Fig. El the flux of the antiproton signal on the top of the atmosphere is depicted 
for various wino masses as a function of the antiproton kinetic energy. In this figure, 
we use the astrophysical parameters of the median set in Table which gives the 
minimal for the B/C analysis |S2|- The parameter of solar modulation is set to 
be 500 MV, which corresponds to almost minimum solar activity. For comparison, 
the background flux is also shown as a dotted line. As shown in this figure, it is 
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Figure 7: Antiproton flux from tlie wino-like neutralino anniliilation and tlie back- 
ground flux on tlie top of atmosphere for three astrophysical parameter sets in Ta- 
ble El Here, we take the wino mass 2 TeV. The solar modulation parameter is set 
to be 500 MV as in the previous figure. 

implausible to exceed the background for almost all region of the wino mass. The 
exception is only the case of the wino mass around 2 TeV. In this case, the wino-like 
neutralinos annihilate resonantly as discussed in section |21 and the signal flux is 
almost comparable to the background flux at high energy (T > 100 GeV). 

Let us discuss the uncertainties in the prediction of the signal antiproton flux. 
In Fig. |7| we show the signal antiproton and the background fluxes on the top of 
atmosphere for three astrophysical parameter sets in Table El in order to see the 
dependence on choice of astrophysical parameter sets. Here, we take the wino mass 
2 TeV and the solar modulation parameter is 500 MV. 

The uncertainties in the astrophysical parameters lead to an uncertainty of a 
factor O(IOO) for the prediction of the signal antiproton flux. The flux from the 
neutralino annihilation depends strongly on the astrophysical parameters, especially, 
the value of the thickness of the diffusion zone L. A larger L means more injection of 
signal antiprotons and leads to larger signal flux. On the other hand, the background 
flux has no strong sensitivity to L, as shown in the figure. As a result, we can observe 
the signal antiprotons only when L is large enough. The situation is very different 
from the positron signal, in which only the DM within a few kpc contributes to the 
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flux due to the rapid energy loss of high-energy electrons. 

Other uncertainties come from the choice of the DM halo profile. The signal an- 
tiprotons in the cosmic rays can travel far from the solar system, while the positrons 
are originated within a few kpc. Therefore, the prediction of the signal antiproton 
flux varies with uncertainties in the halo profile j2HllS2|- This uncertainties depend 
on the size of diffusive halo, especially L. For large 15 kpc) the antiproton 
flux from the neutralino annihilation may be changed by several tens of percent de- 
pending on the choice of halo profile. However, for moderate value of L < 5 kpc, 
the uncertainties of halo profile is negligible with respect to other ones such as the 
diffusive halo size. 

Finally, we concentrate the issue of discrimination of the signal from the back- 
ground. We showed that the wino-like neutralino annihilation with mass around 
2 TeV leads to the signal antiproton flux comparable to or larger than the back- 
ground one. In this bump appears in the antiproton spectrum. However, 
it might be still difficult to recognize presence of the bump from the observed an- 
tiproton spectrum in the upcoming experiments such as PAMELA and AMS-02, 
compared with cases of the positron flux from the DM annihilation. 

The positron flux has a high-energy component, which is produced more directly 
from the leptonic weak boson decays, so that the signal-background ratio for the 
higher energy positron flux is better. However, since the spectrum of the signal 
antiproton flux is featureless even at high energy, it suffers more from uncertainties 
of the background. 

The background antiprotons at high-energy mainly come from the interaction of 
the primary protons in the cosmic rays as discussed in the previous section. The 
proton spectrum is mainly determined by the source term, that is, the injection from 
supernovae. The high-energy protons are produced by the shock-wave acceleration. 
It implies that the spectrum shows the power-law behavior in terms of the energy. 
However, it is difficult to predict the slope of the power law. Thus, the background 
flux also has an ambiguity in the slope at a high-energy range. Varying the slope 
at the source term is immediately reflected to the slope of the antiproton flux at 
the earth. As a result, with lack of the knowledge of the slope of the source term, 
it is difficult to distinguish the signal from the background using only the slope of 
spectrum. 

It may be important to observe the antiprotons with energy around the neutralino 
mass so that the bump may be recognized. When the wino mass is around 300 GeV 
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and the thickness of the diffusion zone L is large, the whole structure of the bump 
might be figured out in the experiments. However, when the mass is around 2 TeV, 
it would be a hard job to detect antiprotons with such a high energy. 

5 Conclusions 

In this article, we have studied indirect detection of the wino-like neutralino DM 
using cosmic ray positron and antiproton observations. Non-perturbative effect en- 
hances the neutralino annihilation cross section when the mass is larger than about 
1 TeV. Especially, when the mass is around 2 TeV, the cross section is enhanced 
significantly due to the resonance effect of the bound state, which is composed of 
the wino-like neutralinos and charginos. In those cases, the cosmic ray positron 
and antiproton fluxes produced by the neutralino annihilation in the galactic halo 
also enhance, and the sensitivities of the upcoming experiments, such as PAMELA 
and AMS-02, are improved for the heavier neutralino DM. It is noticed that the 
relic abundance of the wino-like neutralino in the universe is explained by the ther- 
mal scenario when the mass is around 2 TeV. It might be difficult to study such a 
heavy neutralinos in experiments except for observation of the cosmic rays. Even 
in the direct DM detection, the sensitivity should cover 10~*^''^~^''^cm^ for the spin- 
independent cross section so that the heavy wino-like neutralino is detected 
We have concentrated mainly the heavy wino-like neutralino and have evaluated the 
positron and antiproton fiuxes from the neutralino annihilation using the diffusion 
model. 

We found that both positron and antiproton fluxes increase signiflcantly around 
the resonance (m ~ 2 TeV). However, the positron flux measurement has more 
prospects to detect the heavy wino-like neutralino DM, compared with the antiproton 
one. The signal positron flux exceeds the expected background for the positron 
energy larger than about 100 GeV, and the spectrums in the positron flux and the 
positron fraction are signiflcantly deviated from the background ones. In addition, 
it is plausible that the signal positron spectrum at high energy is less sensitive to 
the astrophysical parameters in the diffusion model or the DM halo proflle, since 
the positrons we observe are produced within a few kpc around the solar system. 
PAMELA and AMS-02 have good sensitivities in a broad region of the positron 
energy 10 GeV E < 270 GeV. Thus, they may distinguish whether the heavy 
wino-like neutralino is the DM. 
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We have also discussed the HEAT anomaly in a positron energy range 4 GeV 
< < 20 GeV. The positron flux from the heavy wino-like neutralino annihilation 
with mass 2 TeV is consistent with it within the experimental error. It is amazing 
that the wino-like neutralino can explain both the DM relic abundance and the 
HEAT anomaly even when the mass is around 2 TeV. 

The antiproton flux from the wino-like neutralino annihilation may be comparable 
to or larger than the expected background for the mass around 2 TeV. However, it 
is strongly dependent on the astrophysical parameters in the diffusion model. In 
addition to it, it might be difficult to discriminate the signal from the background, 
since the antiproton spectrum is featureless. 
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Appendix A: Solution of diffusion equation for positron 
signal 

Here we show how to solve the diffusion equation for the positron signal from the DM 
annihilation. With use of dimensionless parameter e = {E / IGeV), the equation © 
is rewritten under the steady state condition as 



where K{t) = Kq{C + e") and 6(e) = e^/r. The values for Kq, C, a, and r can be 
read off in text. 

We use the cylinder coordinate, so the differentiation is written as = 
+ r~^dr + dl- The source term Q{e,r,z) including the information of the DM 
annihilation is 



where / means the final state of the DM annihilation and {dNe+ / dt) j is the frag- 
mentation function for the final state /. We impose the boundary condition so that 



K{e)V^U (e, r,z) + - {h{e)f,+ (e, r, z)) + g(e, r, z) = , 



(25) 




(26) 
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the density of positron /e+(e, r, z) becomes zero at the surface of the diffusion zone, 
which is given by a cyhnder with radius R and half-height L. 

Due to the boundary condition, it is convenient to expand the density by the 
zeroth-order Bessel function Jq for the coordinate r and by a sine function for z, 



fe+{e,r,z)= ^ An,m{e)Jo(^rj 

m,n=l ^ ^ 



sm(— (^-L)) , (27) 



where Cn are successive zeros of the function Jq. Using the expansion, it is obvious 
that the density satisfies the boundary condition above. 

We comment on the some properties of the Bessel function Jq here. It satisfies a 
following differential equation, 

and has a following orthogonal relation 

£rdrJo (|r) Jo (|r) = , (29) 

where Ji is the first-order Bessel function. 

Substituting Eq. (fTfj) into the diffusion equation and using the differential equa- 
tion and orthogonal relation above, we obtain 

de [r^^ AL^ J "'^ " ^2^n,m[e) . [60) 

Here, we also expand the source term Q by the Bessel and sine functions, and the 
coefficients of the expansion Qn,m are written 

QnAe) = Jfl^^ I dzQ{r, e) Jo (|r) sin [^{z - L)) . (31) 

The boundary condition for An,m is given by An^mi^iaa.^ = 0, where emax = 
max (supp((5n,m)) (~ ^)-^ The condition means that A„ m(e) = if e > emax- The 
function y4„ ^(e) must be a continuous function of e. Solving Eq. (jHUI) . we obtain 



An m / d€ tQyi ryiif- )~^ 



X exp 



Substituting An^m into Eq. ^n\ . we obtain the number density fe+{e,r,z). 
^The symbol, supp(/(x)), implies regions of x in which f{x) ^ 0. 



(32) 
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Appendix B: Solution of diffusion equation for an- 
tiproton signal 



The strategy of solving the diffusion equation for the antiproton signal in Eq. ()15|) 
is essentially the same as that in the positron case. In the equation, the term 
representing the energy loss of the particles, which has differentiation with respect 
to T, is absent. It makes it much easier to solve the equation than that of positrons. 
Thus, we show only the result here. For more detailed calculations, see Refs. P^ I3Uj. 
The number density of antiprotons at the solar system, /p(T, r©), is given by 



oo / 

^exp (- 



y,{z) = 2 / dz'exp { sinh ( ''^ ) g,(T,/) , (33) 



2K 

where Qi is the coefficient of the expansion of the source term Q{T, r, z) by the Bessel 
function, 

oo 

Q(T, r,z) = J2 Q^iT, z)J, {C,-^ . (34) 

i=l 

The parameter Ai in Eq. ()33|1 includes the information about the propagation of 
antiprotons, and it is given as 



A, = 2/iri,ei + K + KSi coth , S, = + . (35) 

From this solution, we can derive the simple relation between the source term and 
the density. Assuming that the source term scales as Q oc E'^", the number density 
behaves as /p(T, r©) oc E^"-^^ , because of A^ oc K{E) oc E^. 
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